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Abstract 

Directed transport of interacting active (self-propelled)Brownian particles is numerically inves¬ 
tigated in confined geometries (entropic barriers). The self-propelled velocity can break thermo¬ 
dynamical equilibrium and induce the directed transport. It is found that the interaction between 
active particles can greatly affect the ratchet transport. For attractive particles, on increasing the 
interaction strength, the average velocity firstly decreases to its minima, then increases, and finally 
decreases to zero. For repulsive particles, when the interaction is very weak, there exists a critical 
interaction at which the average velocity is minimal, nearly tends to zero, however, for the strong 
interaction, the average velocity is independent of the interaction. 
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I. INTRODUCTION 


Diffusion in confined geometries is ubiquitous in nature. The reduction of the coordinates 
in conhned structures can cause the appearance of remarkable entropic effects. More recently, 
physicists have started to study entropic effects in out-of-equilibrium phenomena such as 
transport of particles in corrugated channels l|. Based on the geometry of the channel wall, 
corrugated channels fall into two categories: Compartmentalized channej2-4| and smoothly 


corrugated channels |5l-ll 1| . The relevance of entropic barriers to promote entropic transport 


in conhned environments has been recognized in a variety of situations that include molecular 
transport in zeolites, ionic channels, or in microhuidic devices. The entropic transport in 
these systems yields important and exhibits peculiar properties. 

In previous works, the entropic transport mainly focused on passive particles and few 
works on the entropic transport have involved active Brownian particles. However, active 
matters in biological and physical systems have been studied theoretically and experimen¬ 
tally 12l-l28|. The kinetics of active particles moving in periodic structures could exhibit 


peculiar behaviors. Active particles or agents are assumed to have an internal propulsion 
mechanism, which may use energy from an external source and transform it under non- 
equilibrium conditions into directed motion. Therefore, it would be signihcant important 
to study transport behaviors of active Brownian particles in entropic potentials. Recently, 
Ghosh and coworkers ^studied the transport of self-propelled particles in periodic entropic 
potentials and found that ratcheting current can be orders of magnitude stronger than for 
ordinary thermal potential ratchets. Then, they found that elliptic Janus particles along 


narrow two-dimensional channels can show giant absolute negative mobility 


29 1 and the 


mean exit time of Janus particles in two dimensional cavi 


geometry, particle shape, and self-propulsion strength [30|. Yariv and coworkers [Ml stud¬ 


ies is very sensitive to the cavity 


ied the transport of Brownian swimmers in a periodically corrugated channel by using the 
reduced Fokker-Planck approach. 

In this paper, we mainly studied the ratchet transport of interacting self-propelled par¬ 
ticles in periodic entropic potentials. We focus on hnding how the interaction between 
active particles affects the entropic transport. From numerical simulations, it is found that 
upon variation of the interaction strength, the average velocity exhibits nonmonotonical be¬ 
haviors. For the attractive case, the average velocity hrstly decreases, then increases, and 
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finally decreases to zero. For the repulsive case, the interaction affects the transport only 
for very small interaction strength, for the large repulsive strength, the average velocity is 
independent of the interaction and tends to that in single-particle system. In the regime of 
small interaction strength, there exists a critical value of the repulsive strength at which the 
average velocity takes its minimal value, nearly tends to zero. 

II. MODEL AND METHODS 


In this paper, we consider a set of N interacting self-propelled particles in a periodic two- 
dimensional channel. A self-propelled particle is viewed as characterized by a unit vector 
fii = (cos6*i, sin^i) in the xy plane, dehning the direction of the self-propelled velocity. The 
particles are subjected to both translational and rotational diffusion, with coefficients Dq and 
Do, respectively. The dynamics of the particle i is described by the following overdamped 
Langevin equations. 

rir- «—« f) . -- 

( 1 ) 


^ = vofii - 

* 


dt 




( 2 ) 


where y is the mobility and Vq is the magnitude of the self-propelled velocity. Gaussian white 
noise terms for both the translational and rotational motion are characterized by = 0, 

= Sij6afiS{t - s) and = 0, = SijS{t - s), respectively. Here 

i,j = 1, ...,N label the particles and a, (3 = x,y the coordinates of the space. The symbol 
(...) denotes an ensemble average over the distribution of the random forces. <5 is the Dirac 
delta function. For the sake of simplicity, we have ignored hydrodynamic effects. 

The shape of the channel can be described by its radius u{x) shown in Fig. 1 

where A is the asymmetry parameter of the channel shape and a is the parameter that con¬ 
trols the slope of the channel. The radius at the bottleneck is determined by the parameters 
a, b, and A. 

As for the pair interaction potential V, we consider two cases: (A) the attractive potential 
and (B) the repulsive potential. For case A 


,27rx A 

u(x] = a|sm(—) -|- — sm(^—)J -|- b, 


V{r) = ^kaA, 


( 4 ) 
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FIG. 1: (Color online) Scheme of the entropic ratchet device: interacting self-propelled particles 
moving in a two-dimensional channel. The shape is described by the radius of the channel. 


and for case B 

V{r) = ^, ( 5 ) 

where r is the center to center distance between any two particles, ka and kr are the 
attractive and repnlsive strength, respectively. 

Rectihcation of Brownian particles has been the focus of a concerted effort, both concep¬ 
tual and technological, aimed at establishing net particle transport on a periodic substrate 
in the absence of external biases. The most important quantity characterizing the recti¬ 
hcation in our system is its directional velocity along x direction. Since the Fick-Jacobs 
equation corresponding to the Langevin equations ([I]) and ([2]) can not be solved analytically, 
we have numerically simulated the overdamped two dimensional dynamics of Brownian par¬ 
ticles (fni^ along with the boundary conditions Eq. (|3]) using an improved Euler algorithm. 
In the asymptotic long-time regime, the average velocity of particle i along x direction can 
be obtained from the following formula 



lim 

t—yoo 




( 6 ) 


where 6q is the initial angle of the trajectory. The average velocity after a second average 
over all 9q is 

^ (7) 

The full average velocity is v = b For the convenience of discussion, we dehne the 
scaled average velocity Vs = v/vq through the paper. 


4 









III. RESULTS AND DISCUSSION 


In our simulations, the integration step time At was chosen to be smaller than 10“'^ and 
the total integration time was more than 3 x 10^. The stochastic averages reported above 
were obtained as ensemble averages over 3 x 10^ trajectories with random initial conditions. 
Unless otherwise noted, our simulations are under the parameter sets: « = ^ 

A = 1.0, and N = 4. The simulation results are reported in Figs. 2-6. 
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FIG. 2: (Color online) Average velocity Vg versus the asymmetric parameter A for three cases. 
The other parameters are vq = 5.0, ka = 0.1, ky = 0.001, Dq = 0.1, and Dq = 0.1. 

In Figure 2, we plot the average velocity Vg as a function of the asymmetric parameter 
A for three cases. It is found that the direction of the transport is completely determined 
by the symmetry of the channel. The average velocity is positive for A > 0, zero at A = 0, 
and negative for A < 0. When |A| —>■ 0, the channel is symmetric, so there is no net 
current. When |A| > Ac, the channel is blocked, no particles can pass across the cell of the 
channel. Therefore, there exists an optimal value of |A| at which the average velocity takes 
its maximal value. 

Figure 3(a) shows the average velocity Vg as a function of the attractive strength ka for 
different values of vq. On increasing ka from zero, the average velocity Vg firstly decreases 
to its minimal value, then increases, and hnally decreases to zero. There exist a valley and 
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FIG. 3: (Color online) Average velocity Vg versus the strength of the interaction potential, (a) 
For attractive particles (ka)- (b) For repulsive particles (kr). The other parameters are vq = 5.0, 
Dq = 0.1, and Dq = 0.1 


a peak in the curve and the average velocity takes its minimal value when ka — Vq. This 
features can be explained by the following considerations. The attractive interaction in 
the system can cause two results: (A) reducing the self-propelled driving, which blocks the 
ratchet transport and (B) activating motion in analogy with thermal noise activated motion 
for a single stochastically driven ratchet, which facilitates the ratchet transport. When 
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fca —0, the average velocity attains a constant value (average velocity in noninteracting 
case). When ka < vq, the factor A dominates the transport, so the average velocity Vg 
decreases when ka increases. When ka > Vq, the factor B gradually becomes signihcant 
and the average velocity increases with ka- However, for very large values of ka {ka —)■ cxo) 
all particles are gathered as a single particle, both many-body effects and the individual 
self-propelled driving can be neglected, so the average velocity tends to zero. The average 
velocity of interacting attractive particles is always smaller than that for the noninteracting 
case, which indicates that the attractive interaction always blocks the rectihcation. However, 
for passive Brownian particles, there exist some values of the attractive strength where the 
interaction can facilitate the rectihcation 321. 


Figure 3 (b) describes the average velocity Vg as a function of the repulsive strength 
kr for different values of vq. It is found that on increasing kr, the average velocity Vg 
hrst decreases to nearly zero value, then increases, and hnally tends to a constant. The 
repulsive interaction in the system can also cause two results: (A) reducing the self-propelled 
driving (blocking the ratchet transport) and (B) dispersing Brownian particles (facilitating 
the directed transport). When kr is very small, the factor A dominates the transport and 
the dispersing effect can be neglected. Therefore, the average velocity Vg decreases when 
kr increases from zero. However, when kr becomes large, the dispersing effect becomes 
signihcant, so the average velocity Vg increases. For large values of kr, the dispersing ehect 
completely dominates the transport and the interaction between Brownian particles can be 
ignored. When kr ^ 0 or kr ^ oo, the average velocity Vg tends to the average velocity in 
the single-particle system, therefore, there exits a critical value of kr (very small value) at 
which the average velocity Vg is minimal, and nearly tends to zero. The critical value of kr 
depends on the parameters of the system. 

Figure 4 shows the dependence of the average velocity Vg on the particle number N for 
both attractive and repulsive cases. For the attractive case, the average velocity Vg decreases 
monotonically with increasing N. The average velocity will tend to zero when N ^ oo. The 
is because the effective attractive driving for large number N can be neglected after the 
average . For the repulsive case, the average velocity Vg is independent of the particle 
number N when kr > 0.01. In the inhnitely long channel, the repulsive forces between 
particles disperse particles, the distance between particles become longer and the repulsive 
forces gradually disappear. 
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FIG. 4: (Color online) Average velocity Vg versus particle number N for both attractive and 
repulsive particles. The other parameters are vq = 5.0, D 0 = 0.1, ka = 0.1, = 0.1, and Dq = 0.1 
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FIG. 5: (Color online) Average velocity Vg versus the rotational diffusion Dq for both attractive 
and repulsive particles. The other parameters are vq = 5.0, ka = 0.1, kr = 0.1, and Dq = 0.1. 
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In Fig. 5, we explore the average velocity Vg as a function of the rotational diffusion De for 
both attractive and repulsive cases. When D 0 ^ oo, transport behaviors are similar for two 
cases. The self-propelled velocity changes its direction very fast. The self-propelled velocity 
acts as a zero mean white noise and the nonequilibrium driving in the system disappears, so 
no directed transport occurs and the average velocity tends to zero. For small values of Dq, 
the transport behaviors are different for two cases: (A) For the repulsive case {kr = 0.1), 
the average velocity decreases monotonically when Dq increases. Especially, when D 0 —)■ 0, 
the average velocity tends to a saturate value. In the adiabatic limit, the repulsive force 
can be expressed by two opposite static force Vq and —Vq, yielding the mean zero velocity 
Vg = ['ys(r’o) +'i^s(“'i^o)]/2, which is similar to the single-particle thermal ratchet j33[|. (B) For 
the attractive case, there exists an optimal value of Dq at which the average velocity Vg is 
maximal. The increase of Dq in this case can reduce the self-propelled driving(blocking the 
ratchet transport) and activate the attractive Brownian particles(facilitating the directed 
transport). When D 0 increases from zero, the latter factor dominates the transport and the 
average velocity Vg increases to its maximal value. For further increasing D 0 , the former 
factor takes effect, the average velocity Vg decreases. 




0 . 06 - 

0 . 05 - 

0 . 04 - 

0 . 03 - 


0.02 H 


\ 




V 






0.01 - 


0 . 00 - 


Attractive case k =0.1 

a 

— Repulsive case k^=0.1 


'4 


S 


IE-3 


0.01 


0.1 


D„ 


FIG. 6: (Color online) Average velocity Vg versus the translational diffusion Dq for both attractive 
and repulsive particles. The other parameters are vq = 5.0, ka = kr = 0.1, and D 0 = 0.1. 

Figure 6 illustrates the impact of the translational diffusion Dq on the average velocity Vg 
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for both the attractive and repulsive cases. From the figure, we can see that the curves are 
similar for two cases. When Dq 0, the particles will stay at the bottom of the channel and 
cannot pass the entropic barrier, so the average velocity Vg goes to zero. When Dq ^ oo, the 
translational diffusion is very large, the effect of the asymmetric entropic barrier disappears 
and the average velocity Vg tends to zero. Therefore, there exists an optimal Dq value where 
the average velocity Vg is maximal. 


IV. CONCLUDING REMARKS 

In this paper, we numerically studied the directed transport of interacting self-propelled 
Brownian particles in a two-dimensional periodic channel. It is found that the self-propelled 
velocity acts as the nonequilibrium driving, which can break the themodynamical equilibrium 
and induce the directed transport. The direction of the transport is completely determined 
by the symmetry of the channel shape. The interaction between Brownian particles can 
signihcantly affect the directed transport. For the attractive case: (1) on increasing the 
strength ka from zero, the average velocity Vg first decreases to its minimal value, then 
increases, and hnally decreases to zero; (2) the average velocity Vg decreases monotonically 
with increase of the particle number iV; (3) there exists an optimal value of Dq at which the 
average velocity Vg takes is maximal value, which is different from the noninteracting case, 
where the average velocity decreases monotonically with increase of kr- For the repulsive 
case: (1) the interaction affects the transport only for very small interaction strength, for 
the large repulsive strength, the average velocity is independent of the interaction and tends 
to that in single-particle system. In the regime of small interaction strength, there exists 
a critical value of the repulsive strength at which the average velocity takes its minimal 
value, nearly tends to zero; (2) the average velocity Vg is independent of the particle number 
N when kr > 0.01; (3) the average velocity Vg decreases monotonically with increase of 
the rotational diffusion Dq, especially, it tends to a saturate value when Dq —)■ 0. In 
addition, the average velocity of interacting active Brownian particles is always less than 
that in the noninteracting case (which indicates the interaction always cannot facilitate the 
rectihcation), which is different from the passive case, where the interaction may facilitate 
the rectihcation. 
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